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Using a numerical technique which is different from the iteration method 
of Fox and Li and which is more suitable for the anahjsis of high-order 
modes, we have calculated the diffraction losses and the field distributions 
at tlie reflectors of the low-loss modes of a symmetric confocal resonator for 
Fresnel numbers 0.6 g N m ^ 2.0. We have also computed the modifica- 
tions which result when the two end reflectors are -perturbed by circular out- 
put-coupling apertures centered on the cavity axis. For a range of small 
but useful aperture Fresnel numbers N the aperture diffraction losses can 
be estimated by first-order perturbation theory from the finite-N m results 
appropriate to No = 0. Such estimates fail for those larger Fresnel num- 
bers No for which the mode intensity patterns are significantly distorted 
at the reflectors by the finite coupling apertures. 

I. INTRODUCTION 

Fox and Li 1 demonstrated by numerical iteration that modes in the 
sense of self-reproducing field patterns exist for open Fabry-Perot 
resonators. Using a numerical technique which is different from that of 
Fox and Li and which is more suitable than iteration for the analysis 
of high-order modes, we have calculated the diffraction losses and the 
field distributions at the reflectors of the low-loss modes of a symmetric 
cylindrical confocal resonator for Fresnel numbers 0.6 ^ N m ^ 2.0. 
The results are discussed below. 

An axial section of the symmetric confocal resonator under examina- 
tion is illustrated in Fig. 1. The cavity is bounded at each end by identi- 
cal spherical (parabolic) mirrors whose perfectly reflecting surfaces 
extend over the annular region a ^ p ^ a m . While a comparison of 

333 



334 THE BELL SYSTEM TECHNICAL JOURNAL, FEBRUARY 1965 




Fig. 1 — Axial section of cylindrical confocal laser cavity. The system is sym- 
metric about the cavity midpoint; the two identical reflectors have radii of 
curvature b ; and the reflecting surfaces are confined to the annular region between 
the two radii (a , a m ). 



the a = and ao ^ eigenmodes is instructive as an example of the 
perturbation of eigemnodes by mirror imperfections, this particular 
geometry also has relevance to an aperture output coupling scheme 
proposed by Patel et al. Asymmetric resonators in which, for example, 
only one reflector is pierced by a coupling aperture will be treated in a 
subsequent article. Boyd and Gordon have derived closed-form ex- 
pressions for the eigenvalues and eigenfunctions of symmetric rectangu- 
lar confocal resonators in terms of the angular and radial prolate spheroi- 
dal wave functions. These results were extended to asymmetric 
rectangular confocal systems with output coupling slits by Boyd and 
Kogelnik. 4 Generalized prolate spheroidal functions relevant to the 
cylindrical confocal geometry have been defined by Slepian. Basic 
expressions are summarized in a review article by Kogelnik. 

Assuming that the dimensions of the resonator in Fig. 1 are large 
compared to the wavelength X of light in the cavity, we define the reso- 
nator eigenmodes from the same scalar formulation of Huygens' princi- 
ple used by other authors. 1 ' For the cylindrical confocal geometry the 
field amplitude at the reflectors for a typical mode can be written in 
the form 



Fiv(p,<p) = fip(p) exp (-il<p), 



(1) 



where (p,<p) are radial and angular coordinates in a plane perpendicular 
to the resonator axis and where (l,p) are integral quantum numbers. 
For a symmetric system with identical mirrors, the field amplitude at 
one reflector must be a constant multiple of that at the other. This 
self-reproducing requirement together with Huygens' principle gives 
the following integral equation which must be satisfied by the radial 
function f [p (p ) : 



K lP fiM =^/ ao d P ' P 'Ji(^p)f lJt { P ') 



(2) 
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Here Ji(z) is the Bessel function of order \l\, (a ,a„>) are the p radii 
limiting the reflecting surfaces (Fig. 1), and b is the mirror separation 
and radius of curvature. 

The magiutude of the eigenvalue ki p determines the diffraction loss 
of the (Ip) mode: 

power loss/pass = 1 — | ki p | . (3) 

The phase of the eigenvalue determines the resonant wavelength : 

resonant X = ±irb\ (/ + 1)t - 2 Arg K ip - 2im\~\ (4) 

where n is an arbitrary integer. 

If we normalize the functions fi P (p) over the surface area of the mir- 
rors, then 

% r dp pfiMfM = 8 PQ , (5) 

OX J ao 

where 8 pq is the Kronecker delta symbol (8 pq = 1 if p = q, and d PQ = 
if p y^ q). The orthogonality indicated in (5) for p ^ q follows immedi- 
ately from (2) when the eigenvalues ki p , K iq are nondegenerate and can 
be imposed if they are degenerate. We choose the arbitrary sign of the 
function // p (p) such that f ip (0 + ) > 0. 

For numerical calculations it is useful to replace the radial variable 
p by a dimensionless variable r defined such that 

N(p) =r 2 ^p 2 /\b (0) 

is the Fresnel number appropriate to the radius p. We characterize the 
hole and mirror radii (a , a m ) by Fresnel numbers 



2 

2 an 



»-*•-£■ »--'-• -5- (7) 

In place of the function fi p (p) we introduce a function 

gi P (r) = Zip(r\/X6) (8) 

for which (2) and (5) become: 

Kirfi.ir) = 2tt r m rf/r / J / (27rrr , )o, p (/); (9) 

drrgi p (r)g lq (r). (10) 

■0 

The sign convention /„(0 + ) > requires gi p (0 + ) > 0. 

The eigenvalue equation (9) for the confocal geometry is atypical 
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in the sense that it can be transformed to an equation having a Hermit- 
ian kernel: 2w(rr) i Jt(2irrr). This fact implies that the eigenvalues 
Ki p are real and that the largest eigenvalue (characteristic of the mode 
with lowest loss) can be computed by a variational method. These two 
features do not obtain in other laser geometries for which the integral- 
equation kernel is generally symmetric but not real — that is, not 
Hermitian. 7 While we do therefore expect some qualitative differences 
between the properties of confocal and nonconfocal geometries, we 
can infer from the work of Fox and Li 1 and other authors 8,9 that many 
features are similar. Both the iterative technique of Fox and Li and the 
kernel-expansion-truncation technique we describe below can be ap- 
plied to nonconfocal as well as confocal systems. 

If we assume that the set of functions gi p (r) is complete, we can under- 
stand the iterative method of Fox and Li as follows. Given an arbitrary 
initial field amplitude g w (r), we express it in the form 

g m (r) = ZjJ4„(r). (11) 

Substituting this expression into the right-hand side of (9), we obtain 
on the left-hand side 

{1) (r) = ^pCpKipgipfr), 

the field amplitude after one transit of the optical cavity. Using this 
function on the right-hand side of (9) and repeating this iterative pro- 
cedure, we obtain after n iterations 

g (n) (r) = Z/WV*Cr), (12) 

the field amplitude after n transits. In the limit of large n only terms 
belonging to the eigenvalue of largest magnitude represented (C p ^ 0) 
on the right-hand side of (11) will remain. All other terms will be re- 
duced in proportion to (| ki p \/\ m p | ma x) n - If the two largest eigenvalues 
are sufficiently different, this procedure conveniently yields for each 
angular quantum number I the eigenvalue of largest magnitude, the 
eigenvalue of second-largest magnitude (through the rate of convergence 
of the iteration), and the two amplitude functions belonging to these 
eigenvalues. Results for the cylindrical confocal resonator are given by 
Fox and Li. 1 

In our analysis of (9) we have chosen to apply a different technique 
from that outlined above. Briefly, we expand the Bessel-function kernel 
in (9 ) as a power series 

©l °° I i \m-l / \ 2(m-l) 

S (. + »ii)K— l). G) ' (13) 
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truncate the series after a finite number M of terms, and reduce the 
integral eigenvalue equation (9) to an M -dimensional-matrix eigen- 
value equation which we can easily solve numerically with standard 
matrix-diagonalization routines. [The reduction of (9) to a matrix 
equation is described in the Appendix.] A similar technique has been 
used for the plane cylindrical geometry by She and Heffner. For the 
confocal system for which real -number algebra is sufficient, our com- 
putations require slightly less computer time than the iterative method. 
Moreover, they give the eigenvalues and eigenfunctions of higher-order 
modes, whereas for each I the iterative scheme of Fox and Li is practical 
only for the two largest eigenvalues and their amplitude functions. 

Related methods have been utilized in limited calculations by other 
authors. 9, 10 In place of the power-series expansion (13) it has been 
suggested 10 that one utilize an expansion in associated Laguerre func- 
tions. Because the associated Laguerre functions are the exact eigen- 
functions of the infinite-mirror problem (ao = 0, a m = «>), one might 
expect that fewer terms are required than for the power-series expansion 
(13) and that one can thereby simplify the solution of the matrix 
eigenvalue equation. While such considerations may indeed be relevant 
for Fresnel numbers so large that the matrix eigenvalue problem based 
upon (13) becomes prohibitive, the advantages are largely offset for 
small Fresnel numbers (N m < 4) by the increased effort required to 
compute the necessary overlap integrals. A similar remark applies to 
the Fourier-Bessel expansion used for the plane cylindrical geometry by 
Bergstein and Schachter. 9 

One other numerical technique, different from both the iterative and 
the expansion-truncation techniques, deserves brief mention. If one 
approximates the integral in (9) by a sum over small but finite radial 
intervals, one has in effect reduced the integral equation to a matrix 
eigenvalue equation. If the number of intervals is not too large (<50) 
it is practical to solve this problem directly, although for small Fresnel 
numbers considerably less effort is required with the iteration or ex- 
pansion-truncation techniques. 

In the following section we present results appropriate to the sym- 
metric cylindrical confocal geometry in the absence of coupling apertures 
(AT = 0). We compare those finite-iV m results with expressions derived 
by first-order perturbation theory from the infinite-iV™ eigenfunctions 
and find significant discrepancies. In Section III we indicate how finite 
coupling apertures (No ^ 0) modify these results and derive simple 
mathematical expressions which approximate the machine-computed 
results in useful regions. In Section IV we briefly discuss the far-field 
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output patterns of the aperture-coupled resonator. In the final section 
we briefly recapitulate some of our conclusions. 

II. EIGENVALUES AND EIGENFUNCTIONS WITH NO MIRROR APERTURES 

(No = 0) 

Using the kernal expansion-truncation technique outlined in the 
preceding section and in the Appendix, we have computed the eigen- 
values of (9) for Fresnel numbers N m in the range 0.6 ^ N m ^ 2.0. 
Where they overlap, our results agree with those of Fox and Li 1 and 
other authors. 10 

In our calculation we retained M = max (ION + 1,10) terms of the 
truncated series (13), where, if r ^ r m is the maximum radius of in- 
terest in the field amplitudes gi p (r), we define N = r m r/\b ^ N m . This 
choice insures that the remainder 



Ji(2irrr ) - 2^ 7 r~i TVTT , ,< 

m =i (m + I — l)!(ra — 1)! 



m- If _ M '\l+2(m-l) 

(14) 



will never be greater than 0.001 for the relevant radii. We have indicated 
in Table I for N m = 0.8 and No = how the eigenvalues of the three 
lowest-loss modes converge as the number M of terms increases from 
1 to 10. 

In Fig. 2 we have plotted the power loss/pass = 1 — | ki p | 2 of the 
least lossy modes for 0.6 ^ N m ^ 2.0, No = 0. It is noteworthy that no 
modes have less than 1 per cent loss/pass for N m ^ 0.7, that only two 
modes have such losses for N m ^ 1.0, but that ten modes have less than 
1 per cent loss/pass for N m ^ 2.0. The number of low-loss modes in- 
creases very rapidly for N m > 1.0 so that, whereas N m ^1.0 is in one 
sense a small Fresnel number, N m = 2.0 is already rather large. 

In Figs. 3-8 we have indicated for various low-order modes how the 
field amplitude and intensity varies with radius on the end reflecting 
surfaces for No = and N m = 0.8, 1.6. From the intensity plots it is 
clear that the power loss/pass increases as the mode order increases 
because the higher-order eigenfunctions have more intensity lying out- 
side the reflecting mirrors (and hence lost) than do the low-order eigen- 
functions whose intensity is more concentrated near the mirror center. 

From (9) and (13) it follows [compare (32) in the Appendix] that as 
r -+0 

0ip(r) -> Gi(Zp) [£] r\ (15) 
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Table I 
A. Dependence of Eigenvalues on Number M of Terms 

in Series (13)* 



M 


K00 


K01 


K01 


1 


1.25853309 








2 


0.90311672 


-3.6815805 





3 


1.7198455 


-0.32162843 


0.83714554 


4 


0.98489118 


-1.1885874 


0.24997641 


5 


1.001097(5 


-0.71086465 


0.20030462 


6 


0.99744669 


-0.77432308 


0.12014589 


7 


0.99780308 


-0.76531222 


0.13087265 


8 


0.99777117 


-0.76618958 


0.12962865 


<J 


0.99777343 


-0.76612108 


0.12973537 


10 


0.99777330 


-0.76612543 


0.12972805 



B. Change in Eigenvalues as Number of Terms 
in Series (13) Increases* 



M 


K00 


K01 


KOI 


2-1 


-0.35541637 


-3.6815805 





3-2 


+0.81672878 


+3.35995207 


+0.83714554 


4-3 


-0.73495432 


-0.86695897 


-0.58716913 


5-4 


+0.01620642 


+0.47772275 


-0.04967179 


6-5 


-0.00365091 


-0.06345843 


-0.08015873 


7-6 


+0.00035639 


+0.00901086 


+0.0172076 


8-7 


-0.00003191 


-0.00087736 


-0.00124400 


9-8 


+0.00000226 


+0.00006850 


+0.00010672 


10-9 


-0.00000013 


-0.00000435 


-0.00000732 



* Tables IA unci IB are computed for the case N„> = 0.8, N = 0. 

where Gi(lp) is a constant. This r l dependence of the field amplitude 
and a corresponding /•"' dependence of the intensity is apparent in Figs. 
3-8. Because only the angular-independent (1 = 0) modes have nonzero 
intensity at r = 0, we anticipate that the I = modes are much more 
sensitive to a coupling aperture centered at r = than are the I ^ 
modes, a fact confirmed by the finite-No calculations to be discussed in 
the following section. 

For infinite mirrors without apertures (N m — > °°, Na = 0), 



«lp — ( _ 1)' 



and 



Qipir) = 






(2r7- 2 ) m e— t L p l (2Tr 2 ), 



(16a) 



(16b) 
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Fig. 2 — Power loss/pass versus mirror Fresnel number N m for low-loss modes 
of resonator having no output-coupling aperture (No = 0). 
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Fig. 3 — (a) Field amplitude gi P (r) and (b) field intensity grfir) for modes 
dp) = (00), (01), (02), and (03) with N m = 0.8 and A/ - 0. 
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Fig. 4 — (a) Field amplitude gi P (r) and (b) field intensity gi P 2 (r) for modes 
(/p) = (10), (11), (12), and (13) with N m = 0.8 and N = 0. 
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-1.0 - 




Fig. 5— (ti) Field amplitude gi P (r) and (b) field intensity gi P 3 (r) for modes 
[lp) = (00), (01), (02), and (03) with A',„ = 1.6 and N a = 0. 



. (17) 



where L p ' (z) is the associated Laguerre polynomial 11 

i , , _<i^_ d* , _, P+ i f. (p + !(-«)- 

** w " p! f / 2 p^ e z ' £f (p -m)!(Z + m)!m! 

Low-order Laguerre polynomials are 

U\z) = 1; Li'W -1 + 1 -*; 

I* (*) - *t(I + 2) (I + 1) - 2z(Z + 2) 4- z 2 ]. 

The finite-N,,, results we have computed transform continuously into 
the solutions (10) as N m increases. If for No = and arbitrary iV m the 
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Fig. 6 — (a) Field amplitude gi P (r) and (b) field intensity grfir) for modes 
{lp) = (10), (11), (12), and (13) with N m = 1.6 and N = 0. 

344 



CON FOCAL LASER RESONATOR MODES 



345 




Fig. 7 — (a) Field amplitude (ji,,(r) and (b) field intensity gi P ? (r) for modes 
(lp) = (20), (21), (22), and (23) with JV,„ = 1.6 and iV = 0. 

integer p ^ orders the eigenmodes of a given angular quantum number 
I with respect to increasing power loss/pass, then the eigenvalue 



Kip = (— 1) P | Klj 



(18) 



and the amplitude function has p zeros in the interval < r < r m . 
If we assume for finite N,„ that the low-loss eigenf unctions approxi- 
mate the limiting expressions (16b), we can estimate the deviation of 
k^ from its infinite-A^,„ value (16a) by first-order perturbation theory: 

1 - (-DV (pert) = f dxx l e-*[L p l (x)] 2 . (19) 

J 2TAT m 
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Fig. 8— (a) Field amplitude giJr) and (b) field intensity gi p *(r) for modes 
(lp) = (30), (31), (32), and (33) with N m = 1.6 and JV = 0. 



Two special cases are 
1 — wo (pert) = c~ 2 



1 + koi (pert) = [1 + (2%N m ) 2 ]e' 



We compare these estimates with computed values in Table II. The 
errors between the computed and estimated values in Table II, while 
small relative to the eigenvalues themselves, are nevertheless significant 
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Table II — Deviation of Eigenvalues from Infinite-Mirror 

Values (16b) 



N m f 


1 — KIIO 


1 — >coo(pert)t 


1 - «oo(asym)tt 


1 + K01 


1 + xoi(pert)t 


1 + Koi(asym)tt 


0.6 


1.82 X 10" J 


2.31 X 10-' 


2.52 X 10" 1 


0.545 


0.351 


5.73 


0.8 


2.23 X 10- 3 


0.50 X 10- 3 


2.72 X 10- 3 


0.234 


0.172 


1.099 


1.0 


2.38 X 10~ 4 


1.87 X 10- 3 


2.75 X 10" 4 


0.37 X 10" 2 


7.50 X 10~ 2 


0.174 


1.2 


2.38 X 10-* 


5.31 X 10- 4 


2.08 X 10- 5 


1.20 X 10" 2 


3.07 X 10~ 2 


2.43 X 10- 2 


1.4 


2.24 X 10-" 


1.51 X IQr* 


2.53 X 10~ 6 


1.80 X 10- 3 


1.19 X 10" 2 


3.13 X 10- 3 


1.0 


~2.7 

x 10-7 j 


4.31 X 10" s 


2.34 X 10-' 


2.40 X 10" 4 


4.39 X 10~ s 


3.79 X 10-* 



* N 9 = 0. 

t Estimated values computed from (19) of first-order perturbation theory. 

ft Estimated values computed from asymptotic (20). 

j The accuracy of k 00 for N m = 1.6 is limited by machine rounding errors. 

when compared to the difference 1 — ( — 1) p ki p = 1 — | kj p | which for 
these eigenvalues is roughly one-half the power loss/pass (13). The errors 
arise because the real eigenfunctions are not identical to the limiting 
expressions (16a). Slepian 5 has derived more accurate asymptotic re- 
sults appropriate to N m large : 

r(8TN m )' 



1 - (-l) p /a p (asym) = 

Two special cases are 

1 — koo (asym) = 8ir 2 N m e~ 



t l+2p+l -4irJV„ 



Pl(j + f>)l 



[ 1 + (rJ]- (20) 



1 + Koi(asym) = ir(8wN m )e 



Values computed from these expressions are also listed in Table II. 

In Table III we have listed values at r = of g ip (r) for (Ip) = (00), 
(01), and (02). These values are consistently less than the values 
predicted from the infinite-iV,,, functions (16b) renormalized to the finite 
interval (0,r,„): 

X{1- 



P 



(/ + 



!— f° dx x l e~ x [L p l (x)A 

P)\ J 2rN m j 



(21) 



for which 



g 0p (0) =2 h ll- ^ dx e-*[L p \z))*\ \ (22) 

The differences between the calculated and estimated results again re- 
flect the distortion appropriate to finite N m of the eigenfunctions (16b). 
For a given angular quantum number I, this distortion is generally less 
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Table III — Field Amplitude at Mirror Center for I = Modes 



N m ' 


goo(0) 


goi(0) 


«oi(0) 


0.6 


1.277 


1.225 


1.616 


0.8 


1.321 


1.125 


1.373 


1.0 


1.346 


1.151 


1.193 


1.2 


1.300 


1.209 


1.089 


1.4 


1.369 


1.253 


1.081 


1.6 


1.375 


1.281 


1.130 


2.0 


1.384 


1.315 


1.220 


CO 


-v/2 = 1.414 


1.414 


1.414 



* N = 0. 

significant in low-order modes than it is in high-order modes. (The 
former are usually the only modes relevant to laser oscillators. ) That this 
is so can be understood if we recall that for the cylindrical confocal 
geometry gi p {r) will have p zeros in the interval n < r < r m . As the 
interval (r , r m ) is compressed within the interval (0, °° ) appropriate to 

(16a,b), those functions gi p (r) whose zeros initially lay outside (r , r m ) 
will clearly be more distorted by the compression than will those func- 
tions which have no zeros (p = 0) or those whose zeros lie well within 

(ro,r m ). In Table IV we have tabulated the zeros within (r , r m ) of 
several low-loss eigenfunctions for different N m . 

III. EIGENMODES FOR FINITE MIRROR APERTURES (No ^ 0) 

In this section, in order to distinguish the No = and the No ^ 
results, we mark the eigenvalues and eigenfunctions for No = by a 
superscript "0": Ki p ,gi p (r) a . As in the preceding section, we assign the 
integer p ^ to the No = modes in the order of their increasing 
power loss/pass: | ki p | > | ki pH ° | . We identify the No ^ modes by 



Table IV 


— Zeros of g tp (r) in the 


Domain 


= r <r <r 


'i 






/ = o 





l 


i 


2 


2 


N m * 


r,„t 














/-t = i 


2i 


2. 


l 


2i 


2, 


l 


2i 


2i 


0.6 


0.775 


0.434 


0.318 


0.662 


0.573 


0.432 


0.699 


0.641 


0.502 


0.718 


0.8 


0.894 


0.433 


0.339 


0.735 


0.604 


0.474 


0.792 


0.705 


0.560 


0.821 


1.0 


1.000 


0.425 


0.345 


0.772 


0.606 


0.496 


0.859 


0.734 


0.590 


0.903 


1.2 


1.095 


0.418 


0.342 


0.780 


0.598 


0.504 


0.897 


0.737 


0.621 


0.965 


1.4 


1.183 


0.414 


0.336 


0.774 


0.590 


0.501 


0.909 


0.729 


0.028 


1.002 


1.6 


1.265 


0.412 


0.330 


0.766 


0.585 


0.493 


0.905 


0.721 


0.626 


1.016 


00 


oo 


0.399 


0.305 


0.707 


0.564 


0.449 


0.868 


0.691 


0.564 


0.977 



* iVo = o. 

t The function gi P (r) has p zeroes in r < r < r„ 
% r m = NJ. 
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the indices that those modes would carry if they deformed continuously 
from No = 0. As we shall see, it is not necessary that | ki p \ > | k 1p+ i | 
when JVo^O; however, the No ^ modes do have the properties (15) 
and (18) of the N = modes. In addition, the field amplitude gi P (r) 
will continue to have p zeros in the reflecting interval r < r < r m 
(r = Ofor iVo = 0). 

Using perturbation theory to express the No ^ eigenfunctions in 
terms of the No = functions, we have to first order: 1 



gi P (r) = gi P (r)° {l + tt jH dr' rV P (/)°] 2 } 

/ \0 fro 

' n ici q 2tt / dr r gi p {r ) g lq (r ) . 

o^d K/t. — Kla J 



(23) 



Q^p K/ p Klq 

To second order, the eigenvalues are 



o 

Kip = Kl p 



jl - 2x J ° dr r[gi p (rY 

+ Z' -^^-o W F dr WOWr)']. 

q*p Kl p ° — Klq \_ JO J 



(24) 



The factor multiplying gi p (r)° in (23) is a normalization correction 
compensating for the fact that the gi p (r) are normalized in (10) over 
the interval (r , r m ), whereas the gi P (r)° are normalized over the larger 
interval (0,r m ). The first-order correction to the eigenvalue in the first 
term of (24) decreases the unperturbed eigenvalue ki p by that fraction 
of the unperturbed field intensity which falls on the aperture. 

The second terms in (23) and (24) describe eigenf unction mixing by 
the aperture. The amount of mixing depends upon the eigenvalue 
difference as well as upon the strength of the perturbative coupling. The 
circular apertures, centered on the resonator axis, do not mix modes with 
different angular quantum numbers. Because the signs of the eigenvalues 
alternate as in (18), mode mixing in the symmetric identical-mirror 
cavity is strongest among modes with the same p parity ( — l) p . The 
situation is somewhat different in resonators with dissimilar mirrors 
such as obtains in the apparatus of Patel et al. 2 where only one mirror 
is pierced by the output-coupling aperture. In such systems there is 
significant mixing between even-p and odd-p modes. 

Whereas mode mixing will preclude two eigenvalues from actually 
crossing (if the two modes are coupled by the perturbation), there is, 
because of the sign property (18), no such restriction on the absolute 
values | Kip | and | ki p+ i | or, equivalently, on the diffraction losses of the 
(Ip) and (Ip + 1) modes. For some special values of (No , N m ) one can 
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in fact reverse the power-loss progression of the No = case to give 
I xip+i | > | nip | ; however, it is always true for the identical-mirror 
cylindrical confocal system that 



Kip | > | Klp+2 



(25) 



In more general geometries with more complicated eigenvalue phase 
relations than (18) even this restricted condition can be violated. 

Using the kernel expansion-truncation method outlined in the Appen- 
dix, we have computed the effects of finite coupling apertures on the 
properties of cavity eigenmodes. In Fig. 9 we have indicated for a Fresnel 
number N m = 0.8 how a finite coupling aperture with Fresnel number 
No 9^ affects the loss/pass of the lowest-order modes. In the confocal 
geometry the finite aperture affects only the magnitude of the eigen- 
values; their signs (phases) are still given by (18). In no case do the 
eigenvalues belonging to the same angular quantum number I cross 
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Fig. 9 — Power loss/pass versus aperture Fresnel number JVo for low-loss 
modes with N m = 0.8. Dashed curves (00 E ) and (01 s ) are estimates based on 
(26b). 
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(the circular-aperture perturbations do not couple modes with different 
I), although in Fig. 9 we do see for the (00) and (01) modes a reversal 
of the No = sequence | «oo | > | «oi | in the interval 0.08 < N < 0.16. 

Fig. 9 confirms our previous conjecture that the modes with angular 
quantum number I = are more sensitive to deleterious aperture loss 
than are the I 5* modes. When for N m = 0.8 the area of the aperture 
is only 0.6 per cent of the total mirror area (N /N m = 0.006), the loss/ 
pass of the (00) mode has increased to the point where it equals the 
loss/pass of the (10) mode. For this same hole size the losses of the 
(10) mode are virtually unaffected by the aperture. 

In Fig. 10 we have indicated the intensity distribution of the (00) 
mode for N m = 0.8 and iVo = 0,0.01. Except for the normalization 
correction implicit in the first term of (23), the intensity distribution 
for No = 0.01 is nearly identical to that for No = 0. We conclude for 
N m = 0.8 and iVo < 0.01 that eigenfunction mixing in (23) is unim- 
portant and, as a consequence, that the eigenvalues are accurately given 
by the first-order term of (24). Using the infinite-iV™ functions (16b) 
to approximate gi p (r)° in the first term of (24), we can estimate the 
ratio Kip/ Kip analytically. For I = this estimate can be considerably 
improved if we rcnormalize the infinite-AT m functions (16b) by the 
factor 0o P (O)/\/2 computed from Table III. Doing this, we estimate 



KQp — KOp 



1 " \ WO) ] 2 ] o dx e- x [L p \x)]^ . (26a) 
where ki p ° and gi p (0)° are implicitly dependent upon N m . In the limit 
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Fig. 10 — Field intensity gi P 2 (r) of the mode (Ip) = (00) for N m = 0.8 with 
No = (solid curve) and A r = 0.01 (dashed curve). (N m = r m 2 ; N = r ' i .) 
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of small iVo , for which we can in effect replace gi P (r)° by gi p (Q) in the 
right-hand side of (24), this gives 

k 0p - kop°{1 " TNo[go P (0)°] 2 \. . (26b) 

To the same accuracy, the r l dependence of gi p (r) noted in (15) implies 
that ki v = ki p for I ?± 0. The approximation (26b) has been used to 
compute the dashed curves in Fig. 9. For No < 0.02 the fit to the ma- 
chine-computed curves is excellent. 

The effect of a finite aperture (No 7* 0) on the losses of the low-loss 
modes for N m = 1.6 is indicated in Fig. 11. More modes are shown than 
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in Fig. 9 because at the larger Fresnel number more modes have low 
loss (cf. Fig. 2). Observe for any fixed radial quantum number p that, 
consistent with the r behavior noted in (15), modes with low angular 
quantum number I are more sensitive to a small aperture than are the 
modes with higher angular quantum number. The sequence in Fig. 11 
of upward breaks in the losses of the modes (00), (10), (20), • • • is 
particularly striking, as is that for the modes (01), (11), (21), • • • . 

In Fig. 12 we have redrawn those curves of Fig. 11 which pertain to 
the angular-invariant I = modes. The dashed curves in Fig. 12 derive 
from the approximation (26b), which is here valid only for No < 0.0005 
in the (00) and (02) modes and for N < 0.006 in the (01) mode. 
[For N m = 0.8 it applied to all iVo < 0.02.] The approximation (26b), 
based upon first-order perturbation theory, fails when eigenfunction 
mixing becomes significant. Mixing is strong for the (00) mode when 
the losses of that mode due to the finite aperture approximate the edge 
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losses of the (02) mode — that is, when «oo ~ K02 . If we recall that a 
variational principle applies to the present eigenvalue problem [be- 
cause, in contradistinction to more general resonator kernels, the kernel 
in the integral equation (9) is Hermitian], we can view mode mixing 
as an attempt by the field in the low-loss (00) mode to reduce its in- 
tensity at the aperture and to reduce thereby the total (00) loss. Be- 
cause edge as well as aperture losses contribute to the total loss, the 
deleterious edge losses of the (02) mode preclude appreciable (00) -(02) 
mixing until the aperture losses of the undistorted (00) mode approxi- 
mate the edge losses of the (02) mode. Because edge losses decrease 
rapidly with increasing N m (N m = 1.6 is already quite large), the 
aperture losses required for appreciable mixing decrease rapidly with 
increasing N m . The relevance of mode mixing to the breakdown of 
(26b) is clearly illustrated in Fig. 13, where we show the intensity 
distribution of the (00) mode for N m = 1.6 and N = 0, 0.01, and 0.02. 
[The same aperture Fresnel numbers iVo gave insignificant mode dis- 
tortion for N m = 0.8 (Fig. 10).] 

In Fig. 14 is shown the intensity distribution for three other low-order 
I = modes besides the (00) mode for N m =1.6 and N = 0.01. This 
figure should be compared with Fig. 5b, which shows the intensity 
distribution of the same modes for N„, = 1.6 and No = 0. Note that, 
whereas the intensity at r = of the (00) mode decreased as a result of 
aperture mode mixing, the intensity at r = of the (02) mode increased. 
This increase is reflected in Fig. 12 in the sharp rise of (02) losses as 
the (00) and (02) eigenvalues "repel" for JV > 0.003. 




Fig. 13 — Field intensity g [p 2 (r) of the mode (lp) = (00) for N m = 1.6 with 
N = (solid curve) and iV = 0.01, 0.02 (dashed curves). (N m = rj; N = r 2 .) 
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Fig. 14 — Field intensity gi p 2 (r) of the low-loss I = modes for N m = 1.6 and 
No = 0.01. (N, n = r,„ J ; N = r 2 .) Notice how mode mixing has changed the in- 
tensity distribution near the aperture (0 ^ r < r ) from that in Fig. 5(b) where 
N = 0. 



In Figs. 15, 16, and 17 we have indicated how the power loss/pass of 
the low-loss modes varies when for fixed aperture size No the Fresnel 
number N m changes. Notice in the typical Fig. 15 that the losses of the 
(00) mode decrease as N m increases from No until for N m « 0.7 those 
losses saturate at about 11 per cent/pass, approximately the loss pre- 
dicted from (26b) with 0oo(O)° = s/2 and k 00 ° = 1. As N m increases 
beyond 1.3, mode mixing reduces the losses of the (00) mode as the 
modified intensity distribution avoids both the aperture and the reflector 
edges. While by N m =1.6 the (00) mode again has the lowest loss of 
the I = modes, its total loss is greater than that of certain I ^ modes 
and its intensity distribution (Fig. 13) is considerably different from 
the simple Gaussian of (16b). 

A quantity of interest in the design of lasers with aperture output 
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Fig. 15 — Power loss/pass of low -loss modes versus Fresnel number N„ 
aperture Fresnel number No = 0.01. 



for 



coupling 2 is that value N 0c of No for which the losses of the (00) mode 
equal the losses of the (10) mode. All other things being equal, the 
laser will oscillate in the (00) mode for No < N 0c , whereas for No > 
Noc it will operate in the (10) mode or, for large values of No , in still 
another mode (cf. Fig. 11). In Fig. 18 we have plotted N 0c as a function 
of N m . For N m > 0.6 this curve can be accurately reproduced by the 
following expression based upon (26b ) : 

N 0e = (koo° — Kio°)/7ncoo°[0oo(O) ] 2 - (27) 

This result obtains even for large N m because N 0c decreases so rapidly 
with increasing N m that mode mixing is never relevant. 
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Fig. 16 — Power loss/pass of low -loss modes versus Fresnel number N m for 
aperture Fresnel number No = 0.001. 



IV. FAR-FIELD PATTERNS, APERTURE OUTPUT COUPLING 

If we assume that the useful output coupling of the mode dp) is 
exclusively through one of the small reflector apertures of Fresnel 
number A r , then at a large distance d from the relevant output aperture 
and in a direction making an angle with the cavity axis (see Fig. 19) 
the field amplitude will be proportional to 
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Fig. 17 — Power loss/pass of low-los3 modes versus Fresnel number N m for 
aperture Fresnel number No = 0.0001. 
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(28) 



Here <p is the azimuthal angle used in (1); r is the radial variable defined 
in (6); and gi p (r) is the mirror-field amplitude function (8). The deri- 
vation of (28) from the Fraunhofer formula parallels that of (2) and 
(9). 1 The basic approximation used is 

[d 2 + p — 2pd sin 6 cos ip ]* ftj d — p sin 6 cos <p. (29) 

This approximation is suitable when d y> a ^ p and a /Xd = bNo/d <K 2. 
In the important case for which n is so small (N < N 0c is sufficient) 
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230 OPTICAL 
CAVITY 



Fig. 19 — Geometry appropriate to the analysis of the far-field pattern for 
aperture output coupling of cavity. 

that we can approximate gi p (r) by the lowest-order term (15), we ob- 
tain from (28 ) : 



A lp {e,<p) = l / Qh 



l\ dr l 



-«„ a a Na 



1/2 



a sin 9 \ 



(2irdo sin 6 



)■ 



(30) 



The observed intensity pattern for a pure I mode will be proportional to 



, rid 1 



h.{9#) = 77 



dip 
dr l 



aJVo' /2 ~| 2 2 r; / s -.i T 2 /27ra sine\ .. 

^— cos [l{tp - 8i)]Ji+i ( r ) , (31) 

m o_\ \ a 



r=a d sin 



where 5* is an appropriate phase angle. For I = this gives the cele- 
brated Airy diffraction pattern. 1 



V. SUMMARY REMARKS 

We have computed the eigenmodes of a symmetric cylindrical con- 
focal laser cavity of Fresnel number N m ^2.0 and have determined 
how those modes change when a small circular element centered on the 
axis is removed from each reflector. The calculation methods can easily 
be adapted to cylindrical confocal resonators for which the two mirrors 
and mirror apertures have different sizes. (Results relevant to a coupling 
aperture in only one end reflector will be published in a subsequent 
article.) The basic expansion-truncation methods outlined in Section I 
and in the Appendix are quite general 10 and can usefully be applied to 
nonconfocal geometries for which complex-number computations are 
required when the mirror surfaces are not surfaces of constant phase. 

For the cylindrical confocal geometry the results reported above 
indicate that, while the infinite-JVm functions with appropriate nor- 
malization do approximate the low-order finite- A^ m eigenf unctions, 
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there are significant differences which influence the calculation of both 
aperture and edge diffraction losses. For sufficiently small aperture 
Fresnel numbers No the aperture diffraction losses can be estimated by 
first-order perturbation theory based upon the finite-iV™ eigenf unctions 
(or upon the infinite-N m functions renormalized to the finite-iV m am- 
plitude at r = 0). The value of N for which such first-order calculations 
are valid decreases rapidly as the Fresnel number N m increases, because 
for large apertures the field distributions distort (higher-order per- 
turbation theory) to avoid the aperture. This distortion occurs at 
approximately those values of N and JV m for which an observer at one 
reflector, using light of the relevant wavelength and optics limited by 
the radius r m = N m \ can resolve the aperture of radius r = No at the 
opposite reflector. 
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APPENDIX 

Reduction of Integral Equation (9) to a Matrix Equation 

Truncating the series (13) after M terms and substituting the result 
into (9), we obtain 15 (I = \l\w this Appendix) 



fin M ( 1 \*»— 1 /»,„,' \i+2(m-l) 

K, P g lp (r) = 2* ^ dr r 2, (m + z _ 1)!(w _ 1}! M' > 



; .,„ -/ - i -l(m- 1)! ( 



■ / dr (r ) gi P (r , 

[^J g | f/ ... ( t!Ly"/ u Q-W. 



fei[(7n - l)!(m + /- I)!//!]* 



where 
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2tt 



GMv) - 



[( TO + i-l)!( OT -l)!]» 

(33) 



/: 



j / '/ /2\fn— H-J/2 / '\ 

dr r (irr ) gi p (r ). 

Solving (32) for gi p (r) and substituting this expression into (33), we 
obtain after simple manipulations 

ki p GMv) = E 



ti [(m - l)!(m + I - l)!(fc - l)!(fc + Z - 1)]» 

v [(TrA/,,,)^^- 1 - (TTiVo)' 4 -^ - 1 ] , . 

(I + m + A; — 1) 

We have used the definitions (7) to replace (ro ,r m ) by the Fresnel 
numbers (No , N m ) • 

Equation (34) is a matrix equation which must be solved for the 
eigenvalue ki p and for the M vector components G m (lp) appropriate to 
that eigenvalue. When these latter components are used in (32), one 
obtains the eigenf unction gi p (r) appropriate to the ilf-term truncation 
of (13). The normalization condition (10) on gi p {r) is equivalent to the 
condition 

M 

K lp 8 pg = E (-l) m "U,(Zp)(7 m (Zg) (35) 

on the real vector components G m (lp). The sign condition Re gi p (0 + ) > 
becomes G^lp) ^ where, if Gi(lp) = 0, G 2 (lp) ^ 0, etc. 

In programming the above equations for electronic-computer solu- 
tion, one must insure that at each stage the computations maintain 
sufficient numerical accuracy. The relevance of this remark is clearly 
evident from the fact that, while the Bessel function Ji(z) is of order 
unity for all real z ^ 0, some terms of the series (13) will for z y> 1 
be of order (e/2) 2z /2irz » 1. That is, Ji(z) will be the small difference 
of large numbers and care must be taken to insure that such small 
differences are accurately represented. 

The program utilized to compute the results reported in this paper 
requires a nominal 0.0042 hr. of IBM 7094 running time to compute the 
M different eigenvalues and eigenvectors of (34) for M = 20. Timing 
for other values of M varies roughly as M 3 . 
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